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Abstract 

The diffusion of a system of ferromagnetic dipoles confined in a quasi-one-dimensional parabolic 
trap is studied using Brownian dynamics simulations. We show that the dynamics of the system is 
tunable by an in-plane external homogeneous magnetic field. For a strong applied magnetic field, 
we find that the mobility of the system, the exponent of diffusion and the crossover time between 
different diffusion regimes can be tuned by the orientation of the magnetic field. For weak magnetic 
fields, the exponent of diffusion in the subdiffusive regime is independent of the orientation of the 
external field. 
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I. INTRODUCTION 



The study of magnetic colloids is of great importance both from a theoretical and an ex- 
perimental point of view. Recently there has been an increased interest in the study of the 
structural and the dynamical properties of magnetic confined (in particular on the meso- and 
nano-scale) systems due to the possibility of biomedical [l-3| and engineering applications 
4|. Examples of these magnetic systems are ferrofiuid nanofilms js-?] and magnetorheo- 
logical (MR) fluids. For instance, the translational dynamics of a mesoscopic 3D system 
of permanent magnetic dipoles was studied in Ref. jSj, and it was found that the system 
displays signatures of subdiffusive motion due to the strong suppression of orientational 
fluctuations of the magnetic dipoles by the presence of an homogenous external magnetic 
field. The formation of chains of magnetic dipoles (coagulation effect OMl l| ) is also relevant 
for the dynamical properties of these magnetic systems and may lead to different regimes 
of diffusion. Magnetic clusters of dipolar particles were recently investigated experimentally 
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14l | and they may serve, e.g., as drug delivery mechanisms in biological applications. 



The structural properties of magnetic colloids were recently analyzed experimentally [15 1 
and by means of molecular dynamics simulations [16|, where novel field- induced structural 
transitions were observed in confined ferrofiuid nanofilms. 

In comparison with infinite 3D or 2D systems, confined systems exhibit a particular behav- 
ior due to the competition between the confining potential and the interparticle interaction 
potential. For instance, for a 2D system of repulsive particles confined in a circular parabolic 
potential, previous studies clearly identified the effect of the boundaries on the structural 



and dynamical properties of the system, as well as on the melting |17H22|. Another in- 
teresting possibility of confined systems is realised when the 2D system is subjected to an 
external confining potential (e.g. parabolic) in one direction. The system is called quasi-one- 
dimensional (qlD). Such a qlD system of repulsive interact ing pa rticles self-organize in a 



chain-like structure that was recently studied experimentally 23l-l2q , and through analytical 



and numerical calculations 



27 



3l|. 



Diffusion is strongly modified in confined systems, and may lead to single-file diffusion 
(SFD) 32l440l| . which is directly related to the geometrical constrains imposed by an external 
confining potential. Furthermore, qlD systems can be used as models for the study of collec- 
tive phenomena in low dimensional systems, e.g, vortex matter in type-II superconductors 
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4ll . |42| . colloidal particles 43l. l44j| and dusty plasmas. In addition, the mechanisms o 
transport in narrow channels 



ion 




45l | and DNA manipulation using magnetic particles 
can be studied by modelling qlD systems. 

In the present paper we investigate numerically the properties of a system of ferromagnetic 
dipolar particles confined in a one-dimensional parabolic trap (which models a qlD channel) 
coupled to a thermal bath. The orientation and strength of an in-plane external magnetic 
field B are now control parameters that are able to influence the dynamics of the particles. 
For diluted systems, particles are arranged in a single chain structure in the center of the 
parabolic channel. When B is perpendicular to the channel, the magnetic particles interact 
through a pure repulsive potential. For any other orientation of B, an extra attractive term is 
present in the particle-particle interaction potential. The latter can be dominantly attractive 
or repulsive, depending on the orientation of the external magnetic field. In our numerical 
analysis, we perform extensive Brownian dynamics (BD) simulations and calculate the mean- 
square displacement W{t) of the particles for different parameters which characterizes the 
system. For the case of normal diffusion regime (Einstein or Fickian diffusion), one has 
W{t) = Dot"', where Dq is the "free particle" diffusion coefficient, a is the so-called exponent 
of diffusion (in this case, a = 1.0) and t is time. For values of a ^ 1.0, diffusion is said 
to be anomalous. For instance, in the case of SFD, W{t) = 2Ft°' (with a = 0.5) where 
F is the single-file diffusion mobility factor. We show that the application of an in-plane 
homogeneous external magnetic field leads to different regimes of diffusion depending on the 
orientation and strength of the field. 

This paper is organized as follows. In Sec. II we introduce the model system and the 
numerical methods used in our study, including the approximations and limitations of our 
model. In Sec. Ill, we present and discuss our results. Finally, we present the conclusions 
of our work in Sec. IV. 

II. MODEL AND NUMERICAL METHODS 
A. Model System 

Our system consists of interacting dipolar ferromagnetic particles confined in a quasi- 
one-dimensional (qlD) channel and which is in contact with a thermal bath at absolute 
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temperature T. The pair interaction potential Vpahi^) is given by the sum of the dipole- 
dipole term Vdip(r) and the short-range repulsion Vss{r), such as 
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where Tij is the interparticle separation vector between a pair of particles i and j, is 
the permanent magnetic moment of particle i, a is the diameter of each particle and e is 
an energy parameter which characterizes the short-range repulsion between the particles 
and prevent them from coalescing in a single point 48|. We assume identical particles, i.e., 



/i. The qlD channel is modeled by a parabolic confinement potential defined as 
Vconf = ^^'^yf/'^i where m, u and yi are the mass of each particle, the confinement strength 
(frequency) and the y coordinate of the iih particle, respectively. We also apply an in-plane 
homogenous external magnetic field B, which forms an angle with respect to the x-axis. 
The coupling between the magnetic moment of each particle and the external field is given 
by = /Xj X B. In Fig. [H we show a schematic representation of the system under study 
together with the relevant parameters. 
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FIG. 1: (Color online) Schematic representation of the system. The particles have diameter a 
and dipole moment /Xj, which forms an angle 6i with respect to the x-axis. An in-plane external 
magnetic field B is applied with magnitude B and cj) is the angle between B and the x-axis. 

We assume that the motion of the particles is overdamped which is typical for colloids 
moving in a liquid. The equations of motion for the ith magnetic dipolar particle are 

lii = - 5^[Vi(Kiip + V,,)] - ViKonf + Ut). (2) 

j>i 

= T^ + aUt), (3) 

where rj = Xi'k + yiy is the position vector of particle i and Oi is the angle between the vector 
and the x-axis. Furthermore, 7 is the viscosity of the medium and ^i(t) is a stochastic 
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white-noise with the properties: (i) (^i(t)) = and (ii) {Cim(t)^jn(t')) = 2'ykBTSijSmnS{t — t'), 
where m, n corresponds to the components (x, 6), is the Boltzmann constant and T is 
the absolute temperature of the system. 

Note that the equations of motion ([2]) and ([3]) are not coupled to each other, i.e., we neglect 
ihe coupling between translational and angular velocities of each particle (Magnus force 
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50l|). which is a reasonable approximation for our system. Finally, our model system does 
not take into account hydrodynamic effects (particle- fluid and particle- wall interactions), 
which usually have only a small effect on the qualitative behavior of the diffusion properties, 



as recently demonstrated by Euan-Diaz et al. 5l|]- A similar approach was adopted for a 



purely 2D system in Ref. 52 1. 



B. Numerical Methods 

Before we integrate numerically Eqs. ([2]) and (|3]), it is useful to apply the scaling shown 
in Table [T] in order to write these equations in dimensionless form. In so doing, we introduce 
the unit of time as to = (^"^l/^-, where e = ksTo is the unit of energy and a is the unit of 



length. Moreover, Bq = y^e/a^, fiQ = Vscr^, = (^o)^^ and the dimensionless parameter 
u* = m{ujaY /2e controls the strength of the parabolic confinement potential in the y- 
direction. 

Integrating the dimensionless overdamped equations of motion, we obtain the following 



Ermak-type algorithm 



53| for updating the position (r*) and angle {6*) of particle i during 



the simulation time step At* 



r*(At*) = r*(0) + At*f^^. + /\t*{uj*fg* + V2T*At*^*, 



e*{At*) = e*{0) + At*T*^ + V2T*At*^*, (4) 

where = - ^IK,^ + V*], g* = -V;[iy*y] and r*^ = |K x B*|. Furthermore, V^,^ 
and V^* are given by 

_ ■ Ml 3(K-4)(Mr4) 

*^dip - |3 |5 ) {^) 

K: = 4/|r*/^ (6) 

From this point on we will abandon the * notation and all physical quantities are di- 
mensionless, unless stated otherwise. In our simulations, we use the following parameters: 



TABLE I: Transformation of physical quantities into dimensionless (*) form. 



Physical quantity Variables transformations 

Energy V* = V/e 

Position r* = r j /a 



Time interval At* = At /to 

Nabla V* = aVi 

Frequency uj*=uj/ujq 

Temperature T* = T/Tq 

Magnetic moment //* = Mi / 

Magnetic field B* = B/Bq 

At = 1.0 X lO-*', cu = 10.0, /i = 2.0 and T = 1.0. Note that B and T can be related by the di- 
mensionless parameter rj = which is defined as the ratio between the couphng energy of 
a dipole particle with the external field (/x ■ B) and the thermal energy [ksT). Therefore, by 
setting T = 1.0, the dimensionless strength of the magnetic field corresponds to an "effective 
temperature" control parameter. We also use a simulation box of length (in the x direction) 
Lx = 375.0, and linear density p = N/L^ = 0.8. We choose this value of in order to 
cutoff the interaction potential for distances larger than r = Tc = Lx/2 ^ 187.0, at which 
the interaction energy between a pair of particles is approximately Vdip(r)|r^ ^ 1.0 x 10~^. In 
the X direction, we apply periodic boundary conditions and in the transverse direction, the 
system is confined by the parabolic trap, which is controlled by the parameter u. Note that 
in this work we set a value of u which is large enough to prevent particles from bypassing 



each other, as we demonstrated in a previous study 40[ . This forces the system into a strict 
ID chain of particles in the x direction. The initial configuration of the particles is chosen 
randomly and the system is equilibrated during (1.0 — 5.0) x 10^ simulation time steps. Other 
parameters which characterize the system are the magnitude of the external magnetic field 
(B) and the angle 6 between B and the x-axis. Furthermore, the stochastic white noise 

n 

^j(t) is simulated using the Box-Miiller transformation technique [5^ and in all the results 
presented in this work, the error bars in the plots are smaller than the symbol size. 
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III. INTERACTION POTENTIAL BETWEEN TWO DIPOLES 



Before we study the complete system (i.e., the model described in Sec. [TTl), let us first 
analyze the behavior of the dipole-dipole interaction potential Vdip(r) between two particles 
as a function of (cf. Fig. [1]), assuming that both dipoles are perfectly oriented in the 
direction of the external field. In this case, the interaction potential may be written as 




where 9r (cf. Fig. [T]) is the angle formed between the vector r and the x-axis. We assume 
the simplest case, where 9^ = 0°, which means that particles are forming a perfect one- 
dimensional chain along the x direction. The dependence of Vdip [Eq. ([7])] on the distance 
r between two particles is presented in Fig. [2] for different values of 0. We found that for 
(f) ^ 54°, the interaction potential is dominantly repulsive. On the other hand, for (p ^ 54°, 
the interaction potential has a Lennard- Jones form (e.g, = 0° in Fig. [2]). For small values 
of r, the repulsive term 4|r|~^^ is dominant. For intermediate values of r (1.0 < r < 1.5), 
the particle can be trapped in the potential well due to the presence of the attractive part 
in the interaction potential. For larger distances (r — )■ oo), the interaction vanishes. 




FIG. 2: (Color online) Dipole-dipole interaction potential Vdip(r) [Eq. ([7|)] as a function of the 
distance r between two dipoles and for different values of (j). 
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IV. INFLUENCE OF A STRONG EXTERNAL MAGNETIC FIELD ON DIFFU- 
SION 



The influence of a strong homogeneous external magnetic field on the diffusive properties 
of the model system described in Sec. [TTl will now be investigated. The external field B with 



magnitude B = 100 (which is a typical strong field value used in experiments, see e.g. 55l |) 
forms an angle with respect to the x-axis (cf. Fig. [1]). Note that since we set T = 1.0, the 
parameter rj = 200 > 1, which means thermal fluctuations are weak. We now investigate 
how diffusion depends on 0. 

We will study the diffusive properties of the system through the analysis of the mean- 
square displacement W{t) along the x direction, defined as 

where is the number of particles (we use a typical value of A^=300-900 particles), r is 



an arbitrary time origin 56|, St is the time interval between measurements and {■)r is an 



average over different time origins during the simulation [57l |. 
A. Region (I): 55" < (p < 90° 

First, we consider the external magnetic field perpendicular to the parabolic confinement 
channel, i.e. = 90". In this case, the interaction is purely repulsive [i.e. Vdip(r) oc (1/r)^] 
and the mean-square displacement W(t) [Fig. [3](a)] of the system exhibits a subdiffusive 
regime [single-file diffusion (SFD)], with W{t) = 2FjP'^ for time scales larger than the 
short-time normal diffusion regime (STND), which is characterized by W{t) = Dot \58 \. 
The crossover time tc between these two distinct regimes of diffusion can be estimated 56| 
as the time where the curves Dot and 2Fat^'^ intersect: 

Dot, ^ ^tc^ f ^1 . (9) 



The mean-square displacement for = 90° and = 70° are presented in Figs. [3]^a)-(b), 
respectively. We found that for 55° < < 90°, i.e., when the dipole-dipole interaction is 
purely repulsive (cf. Fig. [2]), W{t) has the following behavior 

, Dnt for t <tc 

w{t) = { \^ (10) 
2F„t°-5 for t > t„ 
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where a straightforward calculation using Eq. (Q gives tc ~ 7.58 x 10~^ {Fa ~ 4.79 x 10~^ 
and Dq ^ 0.110 x 10^^). In this region (I), the crossover time tc and the SFD mobility Fa 
are independent of the value of (p. 



B. Region (II): 0° < c/) < 55° 



For the case of 0° < < 55°, the attractive term present in the dipole-dipole interaction 
potential becomes more relevant with decreasing 0. As a consequence, we expect that the 
diffusion of the dipoles to be affected by the orientation of B. We found that for this region 
(II), the system exhibits the STND followed by a subdiffusive regime, with W{t) = 2Fft(0)t°-^, 
where now tc and Fb depends on the angle and 



W{t) 



Dot for t < tc{(f)) 

2F6((/))tO-6 fort >te(0), 



with tc{(f)) ~ {2Fb{(j)) / Dq)'^'^ . The mean-square displacement for 
presented in Figs. [3](c)-(d), respectively. 
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FIG. 3: (Color online) Log-log plot of the mean-square displacement (solid black curves) W{t) as 
a function of the time t ior B = 100 and (a) = 90°, (b) (j) = 70°, (c) (j) = 50° and (d) = 0°. 
The dashed orange lines are a guide to the eye and the crossover time tc for each case is indicated 
by the horizontal arrow. 
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In Figs. Ill^a)-(b) we show the mobihty in region (II) and the crossover time tc as a 

function of 0, respectively. Note that both Fb and tc decreases with decreasing in region 
(II). On the other hand, as stated above, the crossover time tc is constant in region (I). 
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FIG. 4: (a) Mobility Fi, in region (II) as a function of (p and (b) crossover time tc between the 
STND regime and the subdiffusive regime as a function of (j). The solid lines are a guide to the eye. 



The decrease of tc and Fj,, in region (II), with decreasing (p can be explained by the 
decrease of the minimum interparticle distance between neighbor particles [cf. Fig. |5]^a)]. 
When the interaction potential is dominated by the repulsive part of the potential [i.e., region 
(I)], the particles are distributed homogenously along the unconfined direction [Fig. |6]^a)], 
i.e., the minimum interparticle distance between neighbors is approximately constant. In 
region (II), the attractive term in the interaction potential becomes more relevant, and the 
system starts to form clusters of chains. Therefore, the particles are no longer homogenously 
distributed along the channel. The minimum interparticle distance decreases with decreasing 
and the crossover time tc is smaller than in region (I) because particles "feel" the interaction 
with neighboring particles much faster. Also, since the particles can be trapped inside the 
clusters of chains, the mobility (Fh) is reduced with decreasing 0. 

C. Exponent of diffusion (a) in the intermediate (ITR) subdiffusive regime 

In the previous section we showed that the MSD [W^(t)] exhibits two different regimes of 
subdiffusion depending on the region [(I) or (II)]: the exponent of diffusion (a) in the sub- 
diffusive regime changes from a = 0.5 to a = 0.6 as the angle is decreased from = 90" 
to = 0". The exponent a is calculated by fitting the MSD of our simulation data in the 
region of interest (for instance, the ITR regime) according to the relation W{t) oc The 
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FIG. 5: (a) Minimum interparticle distance d between neighboring particles for B = 100 and 
T = 1.0 as a function of the orientation (p of the external field, (b) Exponent of diffusion (a) as a 
function of the orientation (j) of the external magnetic field. Note that d decreases with decreasing 
(f) in the region 0° < (/> < 55°, which is the same region where we found the increase of the diffusion 
mechanism [cf. panel (b)]. The solid hnes are a guide to the eye. 

increase in the diffusion mechanism can be seen in Fig. [5](b), where a is presented as a func- 
tion of the orientation 0. Note that a increases with decreasing 0, which can be understood 
in terms of the dipole-dipole interaction dependence on (p. For (p > 55°, the interaction 
potential is mainly repulsive and therefore it leads the system into a subdiffusive behavior, 
where a = 0.5. The scaling W{t) oc t^'^ has been observed experimentally in repulsive in- 



;eracting particles 35j, and was also found from simulations 4J, |56| and through analytical 
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62| calculations. In this case, the minimum interparticle distance is approximately equal 
to d ^ {p)~^ ~ 1.2. On the other hand, for (p < 55°, the interaction potential exhibits a 
competition between a repulsive and an attractive term (cf . Fig. [2]) . The attractive part of 
the potential forces the formation of clusters of chains [Fig. El^b)], resulting in empty spaces 
along the unconfined direction. This is illustrated in Fig. |5]^a), where the minimum distance 
between particles d is shown as a function of (p>. Note that d decreases with decreasing (p. 
Since the system has a fixed density p, the empty spaces between the clusters of chains 
results in an increase of diffusion, which subsequently gives an exponent of diffusion a that 
is slightly larger than 0.5. 
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FIG. 6: Typical snapshots of the system after 10^ simulation time steps for (a) (p = 90° and (b) 
(J) = 30°. Other parameters are B = 100 and T = 1.0. 

In order to better understand the increase of the exponent of diffusion a, we calculate 
the mean-square displacement of each jth particle [Vrj(t)] using an expression similar to Eq. 

W,it) = {[x,iT + 5t)-x,iT)]')^, (12) 

where j = 1,...,N represents each individual particle and (■),- is an average over different 
time origins during the simulation. In Figs. [71^a)-(b) we show W{t) (open black circles) and 
Wj{t) (gray triangles) for = 90° and = 0°. Note that for the case = 90°, Wj{t) deviates 
very little from the mean-square displacement of the system W{t). In this case the particles 
in the system are distributed homogenously along the unconfined direction. Therefore, the 
diffusion of a tagged particle should be the same as the diffusion of the whole system. On 
the other hand, for the case of = 0°, Wj{t) deviates [much more] from W{t) [than in the 
case = 90°]. This is caused by the asymmetry along the unconfined direction. In this case, 
it is possible that a tagged particle can diffuse differently than the whole system because of 
the formation of clusters of chains [cf. Fig. [6](b)]. For instance, particles which are located at 
the borders of the cluster of chains diffuse faster than particles which are inside the cluster. 
This is the reason for an exponent a that is slightly larger than 0.5 in the case where the 
interaction potential has both repulsive and attractive terms [Region (II), see Sec. IIV Bj . 
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FIG. 7: (Color online) Mean-square displacement of the system [open black circles, and 
mean-square displacement of individual particles [gray triangles, VFj(t)] as a function of the time t 
for two different values of (p = (a) 90° and (b) 0°. The dashed orange lines are a guide to the eye. 
Other parameters are B = 100 and T = 1.0. 

V. WEAK MAGNETIC FIELDS 

In the previous section, we showed that the diffusion mechanism of the system is affected 
by the orientation of the strong external magnetic field. Now, we turn to the question of 
how the magnitude of B influences the diffusive properties of the system. To this end, we 
perform similar simulations using the same parameters of the previous section, but with a 
weaker magnetic field B = 0.1. Note that since we set T = 1.0, the parameter t] = 0.2 < 1, 
which means thermal fluctuations are strong. The mean-square displacement (in log-log 
scale) as a function of the time is presented in Fig. [8] for different values of 0. 
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FIG. 8: (Color online) Log-log plot of the mean-square displacement (solid black curves) W{t) as 
a function of the time t for S = 0.1 and (a) </) = 90°, (b) </> = 70°, (c) c/) = 50° and (d) = 0°. The 
dashed orange lines are a guide to the eye and the crossover time tc for each case is indicated by 
the horizontal arrow. 



There are two important observations regarding the results for 5 = 0.1: (i) note that the 
ITR regime for this case is shifted to larger time intervals as compared to the previous case 
(see Fig. [3]), which is a consequence of the weaker coupling of the dipoles with the external 
magnetic field, leading the system to larger relaxation (crossover) times. Here, the ITR 
regime can be identified in the time interval 10*^ < t < 10^; (ii) since the external magnetic 
field is small (compared to the case of the previous section, B = 100), the coupling between 
the magnetic dipoles and the external field is weaker which results in an approximately (p- 
independent regime of diffusion [cf. Fig. |8]. This means that the exponent of diffusion a 
in the ITR regime is a constant {a = 0.35) which is independent of the orientation of the 
external magnetic field. We will further discuss this particular value of a in the following 
section. Note that, as opposed to the case of strong magnetic field, the clustering of particles 
in a chain-like configuration along the unconfined direction is less pronounced, as illustrated 
in Fig. [9l Note that the orientation of the dipoles of the ferromagnetic particles is almost 
random. 
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FIG. 9: Typical snapshots of the system after 10^ simulation time steps for (a) (/> = 90° and (b) 
(j) = 30°. Other parameters are B = Q.\ and T = 1.0. 

VI. INFLUENCE OF THE STRENGTH OF THE MAGNETIC FIELD 

In this section, we further investigate how the strength B of the external magnetic field 
influences the diffusion of the system. We analyze the case for = 90°, where the SFD 
is found in the ITR regime for B = 100 [see Fig. [5](b)]. From the calculations of the 
MSD using Eq. (|8]) for different values of B, we found that for B > 1.5, the SFD regime 
is always present in the ITR regime, i.e., W{t) oc t^'^. Therefore, we only investigate the 
region 0.1 < -B < 2, and the results are plotted in Figs. [Tor a)-(d). For B = 1.5 [Fig. [TOT a)]. 
as stated above, the SFD regime is present in the ITR regime, which means a = 0.5. 

We found that by decreasing the value of B, the exponent of diffusion (a) decreases from 
a = 0.5 to a = 0.35, as shown in Fig. [TTl The reason for this behavior is explained by 
the following: as the magnetic field is decreased, the coupling of it with the dipoles also 
decreases, leading to an increase in the rotational movement of the dipoles. Therefore, the 
energy of a dipole is distributed between translational and rotational motion. Recall that 
for large values of B (= 100), the dipoles were almost completely aligned with the field. 
The increase in the rotation of the dipoles thus leads to a slowing down of the translational 
diffusion, i.e., a decreases with decreasing B. 

In order to strengthen this conclusion, we calculate the mean-square angular displacement 
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FIG. 10: (Color online) Log-log plot of the mean-square displacement (solid black curves) W{t) as 
a function of the time t for (f) = 90° and (a) B = 2, (h) B = 0.8, (c) B = 0.5 and (d) 5 = 0.1. The 
dashed orange lines are a guide to the eye. 




FIG. 11: Exponent of diffusion a as a function of the strength B of the external magnetic field for 
(j) = 90". The solid line is a guide to the eye. 



(MSAD) Wrotit), which is defined similary to Eq. (E]): 

= ^t^'^^ + ~ ^'^^^^') ' ^^^^ 

where di is the angular coordinate of the ith particle (cf. Fig[T]). The results of calculations 
of the MSAD are shown in Fig. [12] for different values of B. Note that for all values of 
B, Wj-ot{t) saturates after the initial motion. Furthermore, the MSAD curves increases 
with decreasing B, which indicates that the rotational motion of the dipoles increases with 
decreasing strength of the external magnetic field. 



16 



-4 -3 -2 -1 1 

10 10 10 10 10 10 

t 

FIG. 12: (Color online) Log-log plot of the mean-square angular displacement Wrot{t) as a function 
of the time t for (p = 90° and B = 2, B = 0.8, B = 0.5 and B = 0.1. The dotted orange horizontal 
lines correpond to the saturation values of Wrot{t). 

VII. CONCLUSIONS 

We studied a system of interacting ferromagnetic dipoles, confined in a qlD channel, 
that are subjected to a homogeneous external magnetic field. The analysis of the mean- 
square displacement W(t) indicates that the diffusive properties of the system depends on 
the orientation and on the strength of the external field. For the case of strong magnetic 
fields (we considered B = 100 as an example), we found that the exponent of diffusion a 
increases with decreasing orientation [cf. Fig. |5]^b)] of the external field, i.e., directing the 
magnetic field towards the direction parallel to the channel. This increase of diffusion was 
explained by the dependence of the dipole-dipole interaction potential on 0. For > 55°, 
the interaction is dominantly repulsive, leading the system into subdiffusive motion in the 
ITR regime. On the other hand, for < 55°, the interaction potential has a Lennard- Jones 
form, which creates a competition between the repulsive and the attractive term of the 
dipole-dipole potential. The attractive part of the potential leads the system into clusters 
of chains [Fig. |6|b)]. The empty spaces in the system allow for an increase in diffusion. 

For small values of the magnetic field (e.g. B = 0.1), the coupling between the magnetic 
dipoles and B is weak and the dynamic behavior of the system becomes almost independent 
of the orientation of B. This results in an exponent a, in the subdiffusive regime, that is a 
constant [a = 0.35) for all values of the orientation of the external magnetic field. The fact 
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that for weak magnetic fields the exponent of diffusion is smaller than 0.5 (i.e. the slowing 
down of translational diffusion) was explained by the weak coupling of the dipoles with the 
external field, leading to an increase in the rotational motion of the dipoles. 

Our results show that the diffusion mechanism in this system can be controlled by tuning 
the orientation and the strength of the external magnetic field. This will allow one to control 
the dynamics of magnetic particles in narrow channels by simply tuning the parameters which 
regulate the external magnetic field. 
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